Astronomy & Astrophysics manuscript no. Recio-Blanco 


© ESO 2008 


February 5, 2008 





MULTIVARIATE ANALYSIS OF GLOBULAR CLUSTER'S 
. . HORIZONTAL BRANCH MORPHOLOGY: searching for the 
o ■ second parameter. * 

<N : 

^ . A. Recio-Blanco^'^, A. Aparicio^'^, G. Piotto^, F. De Angeli^'^, S.G. Djorgovski^ 

o : 

^ ■ ^ Dpt. Cassiopee, UMR 6202, Observatoire de la Cote d'Azur, B.P. 4229, 06304 Nice Cedex 04, France 
I , e-mail: arecio@obs-nice.fr 

' ^ Departamento de Astrofisica, Universidad de La Laguna, Via Lactea s/n, 38200, La Laguna, Tenerife, Spain 
e-mail: aiitapaj@iac.es 

, ^ Dipartimento di Astronomia, Universita di Padova, Vicolo dell'Osservatorio 2, 1-35122 Padova, Italy 
^ ' e-mail: recio,piotto,deangeli@pd. astro, it 

' ■* Astronomy Department, California Institute of Technology, MC 105-24, Pasadena, CA 91125. 
^ 1 ^ . e-mail: george@astro.caltech.edu 

Institute of Astronomy, Madingley Rd, CB3 OHA Cambridge, UK 
' , e-mail: fda@ast.cain.ac.uk 
• Instituto de Astrofisica de Canarias, Via Lactea s/n, 38200 La Laguna, Tenerife, Spain 

o '■ 

^ ■ Received 8th March 2005; accepted 16th November 2005 
Q ■ ABSTRACT 

C/j Aims. The interpretation of globular cluster horizontal branch (HB) morphology is a classical problem that can significantly 

blur our understanding of stellar populations. 
^ Methods. In this paper, we present a new multivariate analysis connecting the effective temperature extent of the HB with other 
• i-H , cluster parameters. The work is based on Hubble Space Telescope photometry of 54 Galactic globular clusters. 

' Results. The present study reveals an important role of the total mass of the globular cluster on its HB morphology. More 
I massive clusters tend to have HBs more extended to higher temperatures. For a set of three input variables including the 
. . . temperature extension of the HB, [Fe/H] and Mv, the first two eigenvectors account for the 90% of the total sample variance. 
Conclusions. Possible effects of cluster self-pollution on HB morphology, eventually stronger in more massive clusters, could 
explain the results here derived. 
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1. Introduction 

Globular clusters (GC), comprised of chemically homo- 
geneous and coeval populations of stars, represent ex- 
cellent systems for testing stellar models. The various 
sequences that appear in the colour-magnitude diagram 
(CMD) of a globular cluster can be compared to the pre- 
dicted isocrones and theoretical loci. In this way, the prop- 
erties of stars at different stages of evolution, and the fun- 
damental characteristics of the clusters themselves, such 
as cluster distance and age, can be derived. Hence, it is 
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not surprising that the study of GCs has played a pivotal 
role in the development of stellar astrophysics. 

In this paper, we will focus on one evolutionary stage: 
the horizontal branch. HB stars are characterized by core- 
helium burning and shell-hydrogen burning. The star 
location in effective temperature along the Zero Age 
Horizontal Branch (ZAHB) depends on virtually all stel- 
lar parameters (composition, age, rotation, etc., see e.g. 
Rood 1973). The wide colour distribution of the HB, called 
the HB morphology, is the result of large differences in the 
envelope mass of stars having the same core mass, at the 
same evolutionary stage. The HB phase behaves as an am- 
plifier, displaying the record of both initial conditions and 
of any variations and perturbations in the evolution of the 
star from its birth up to the HB stage. Therefore, reading 
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properly the HB morphologies can yield a better under- 
standing of Population II stellar evolution in general, and 
of the specific stellar systems, clusters or galaxies, in par- 
ticular. 

However, it appears that our comprehension of the 
HB phase and its precursors is incomplete. Canonical stel- 
lar theory can not adequately explain the wide variety of 
HB morphologies observed in Galactic GCs, ranging from 
short red HBs to long extended "blue tails" . In particular, 
blue tails still represent a puzzle in the stellar evolution 
model, in the sense that we know what the stars in blue 
tails are, but we do not know how stellar evolution can 
create them. 

To a first approximation, the different temperature ex- 
tension and morphology of the observed HBs were inter- 
preted in terms of the metal abundance variation, the first 
parameter, metal-rich clusters tend to have short red HBs, 
while metal-poor ones exhibit predominantly blue HBs. 
Nevertheless, the previous approximation turned out to 
be too rough very soon. Some other parameter (or set of 
parameters) was evidently also at work, as clusters with 
nearly identical metallicities could show very different HB 
colour distributions (van den Berg 1967, and Sandagc & 
Wildey 1967). One classical example is the pair formed by 
M3 and M13, with [Fe/H] = -1.57 and [Fe/H] = -1.54 
(Harris 1996, in its revised version of 2003) respectively, 
but very different HB morphologies. The variety of pro- 
posed candidates range from cluster age, to helium mix- 
ing, [CNO/Fc] abmidancc, cluster concentration, stellar 
rotation, planets... In fact, the second parameter problem 
has already been the object of an extensive list of studies 
that is important to acknowledge. An increasing amount 
of observational data progressively revealed the complex- 
ity of the scenario (Kraft 1979, Freeman & Norris 1981). 
Important work was developed by the Yale group (e.g. Lee 
et al. 1987, 1988, 1990, Sarajedini & King 1989) interpret- 
ing cluster age as a global second parameter, in the frame- 
work of a Galactic formation picture a la Searlc & Zinn 
(1978). On the other hand, relevant questions outside this 
scenario were also discussed by different authors (Renzini 
& Fusi Pccci 1988, Rood & Crocker 1989, Buonanno et 
al. 1989, Fusi Pecci et al. 1990, Fusi Pecci et al. 1993). 
In particular, the idea that production of hot HB stars 
may be somehow influenced by the dynamical processes 
in the cluster was also carefully explored. In agreement 
with this, Fusi Pecci et al. (1993) and Buonanno et al. 
(1997) suggested environemcnt density as a possible sec- 
ond paramenter. Finally, among other more recent works 
those of Soker & Harpaz (2000) and Catelan et al. (2001) 
can be cited. 

In this work, we have analyzed a homogeneous 
database of 54 globular cluster CMDs from Hubble Space 
Telescope photometry in an attempt to quantify the differ- 
ent dependences of HB morphology on cluster parameters: 
metallicity, concentration, distance to the Galactic center, 
total mass, etc., in other words, to search for the so called 
HB second parameter (s). The data, reduced and treated 
in an uniform way, represent an exceptional opportunity, 



from the statistical point of view, to investigate how HB 
morphology depends on globular cluster properties. On 
the other hand, the multidimensional data set of Galactic 
globular clusters spans a large range in many of their prop- 
erties such as luminosity, metallicity, etc... Therefore, in 
order to reveal the possible complex correlations hidden 
in the HB second parameter problem, we decided to ap- 
ply a multivariate statistical analysis. This approach can 
not only confirm or reveal new correlations, but offers the 
possibility to estimate the relative importance of the vari- 
ous HB dependences and the degree of explanation of HB 
morphology that can be obtained by their combination. A 
similarly motivated study was already done by Fusi Pecci 
et al. (1993). We believe that our new, HST-based data 
set warrants a fresh look at the problem. 

2. The database 

The analysis presented here is based on a Hubble Space 
Telescope snapshot program aimed at mapping the cores 
of all GCs with (m-M)^ < 18.0, using the Wide Field/ 
Planetary Camera 2 (WFPC2) of HST. All the photo- 
metric data come from HST/WFPC2 observations in the 
i^439W and i^555W bands, the WFPC2 equivalents of the 
B and V filters, which are suited for a generic survey and 
constitute the best choice to identify new anomalous HBs. 
In all cases, the PC camera was centered on the cluster 
center. 

The colour-magnitude diagrams and the photometry 
derived from that program have been already published by 
Piotto et al. (2002). Moreover, the database has already 
given rise to a number of works attacking still-open topics 
on evolved stars in GCs (Piotto et al. 1999, Palmieri et 
al. 2002, Raimondo et al. 2002, ZoccaH et al. 1999, Bono 
et al. 2001, ZoccaH et al. 2000, Cassisi et al. 2001, Riello 
et al. 2003, Recio-Blanco et. al. 2004, Piotto et al. 2004, 
Salaris et al. 2004, De Angeli et al. 2005 ). The complete 
database consists in a total of 74 GCs (53 snapshot plus 
21 archive data). For this work, only those clusters whose 
CMD had a well populated HB and enough photometric 
precision to offer reliable estimations of HB temperatures, 
54 of them, were used. 

Due to the severe stellar crowding, the excellent resolv- 
ing power of HST is crucial for the photometrical stud- 
ies of the globular cluster central regions, where accurate 
ground based observations are precluded. Moreover, for 
GCs near the Galactic center, only the central regions are 
sufficiently uncontaminated by field stars to allow a good 
study of the colour-magnitude diagram. 

3. Considered cluster parameters 

3.1. The morphology parameter: maximum effective 
temperature along the HB 

The first step in this work consists in the evaluation of 
the HB morphology of each cluster. In order to have a 
quantitative measure of their extension, we determined 
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the highest effective temperature reached by the stars in 
the HBs of all the clusters in our data set by fitting a Zero 
Age Horizontal Branch (ZAHB) model to the observed 
CMDs. In this way, we can study how the extension of 
the horizontal branch varies with cluster parameters, ft is 
worth to mentioning, however, that HB bimodality, i.e. the 
presence of both a red HB and a blue tail as in NGC2808, 
is not been taken into account with this approach. 

ZAHB models from Cassisi et al. (1999) were fitted to 
our F439W^, F555W CMDs using the values of distance 
modulus and reddening derived in our previous paper, 
Recio-Blanco et al. 2004, for each of the clusters in the 
data set. 

This procedure allowed us to evaluate the highest T^g 

reached by the globular cluster HB and therefore its tem- 
perature extension, as illustrated in Figure 1 for the case 
of NGC 1904 for which the corresponding temperatures 
along the HB arc marked. The errors in this temperature 
determination are difficult to estimate as they depend not 
only on the errors in the distance modulus, {m—M)F555w , 
and reddening, but also on the number of stars in the HB 
and the temperature range we have to deal with. As a 
consequence, the largest errors occur for the smallest low 
central concentrated clusters, and for the most extended 
HBs, where the large bolo metric correction in these pho- 
tometric bands precludes an accurate estimation of T^g. 
However, although the errors can be rather important, the 
general trend of HB morphology with cluster parameters 
is not dramatically affected, as we will see in the next 
Section. 

Finally, we have performed a comparison between our 
Teff parameter, in its logarithmic form, and the Lj param- 
eter from Fusi Pccci ct al. (1993). As explained in their 
Section 3.3.1, the parameter measures the total lenght 
of the HB from the HB red endpoint. The result of the 
comparison is presented in Figure 2. As indicated, the de- 
rived correlation is 0.65. In fact, although there is a clear 
common trend, the spread of the points is rather high, 
most probably due to the difference in the available pho- 
tometric sources. Fusi Pecci et al. (1993) measurements 
come from an extensive, although not homogeneous, col- 
lection of ground based CMDs from 1965 to 1992. On the 
contrary, our measurements were performed using recent 
homogeneous HST data base. Therefore, alghouth both 
parameters rely on a somewhat subjective estimation of 
the terminal HB point, an important part of the observed 
rms is probably coming from the photometric data. 

3.2. Other parameters 

In order to disentangle the dependence of the HB mor- 
phology on as many cluster parameters as possible, we 
have considered the 15 quantities listed in Table 1 : max- 
imum effective temperature along the HB (log(TeffiirB)); 
cluster metallicity {[Fe/H]); total luminosity (My); col- 
lisional parameter (Tcoi)', logarithm of central luminosity 
density in Solar luminosities per cubic parsec (po); cen- 



tral concentration (c=log(rt/rc)); distance from Galactic 
center in kpc assuming Ro=8.0 kpc (Rg-c); Galactic longi- 
tude (1) and latitude (b); core radius (rc); half-light radius 
(r/j); the logarithm of core relaxation time (t^c); the log- 
arithm of relaxation time at half-light radius (trh)', the 
central surface brightness (nv) and the age in a relative 
scale (AgCrei). 

The collisional parameter, listed in column 5, is defined 
as the probability of collisions, per unit time, for one star 
in the cluster, and it was derivated via the formula (King, 
2002): 

Tool = (l0g[5 • W-^^V^^c])/N,tar 

where ic is the core radius in units of parsecs and a is 
the central surface brightness in units of Lq/pc^: 

^ = lQ[-0.i-{l^v-26Al)] 

The total number of stars in the cluster, l^star has 
been estimated by using the integrated visual absolute 
magnitude of the cluster, assuming M/L = 2 and a typical 
mass for the colliding stars of 0.4 mQ. 

Column 16 gives the cluster relative ages from a sub- 
sample of 47 clusters in common with De Angeli et al. 
(2005). They used the so-called vertical method to esti- 
mate ages for a sub-sample of 41 cluster from the snapshot 
database and 30 clusters from the ground based database 
presented in Rosenberg et al. (2000a, 2000b). Sixteen clus- 
ters were in common among the two database and were 
used to assess the consistency of the two catalogs. Our 
analysis will include 3 clusters from the ground based cat- 
alog and 39 clusters from the snapshot one. Five more 
snapshot clusters had age estimates that did not match 
the accuracy of the other determinations and for this rea- 
son were not included in the final version of De Angeli 
et al. (2005), although their ages had been determined 
homogeneously with respect to the rest of the catalog. 
Nevertheless, we decided to include them in our analysis 
given their statistical value. 

The other quantities in Table 1 (see table's caption) are 
taken from Harris (1996, in its revised version of 2003). 

4. Monovariate correlations 

We first explore the simple, pairwise correlations between 
the HB extension and a number of selected globular clus- 
ter properties, via the Pearson coefficient, r. This coef- 
ficient gives the ratio between the observed covariance 
and the maximum possible positive covariance for the two 
evaluated quantities, x and y . Therefore, the value of r 
goes from perfect negative correlation (r= —1) to perfect 
positive correlation (r= -|-1). The midpoint of this range, 
r= 0.0, corresponds to a complete absence of correlation. 

The coefficient of determination, r^, represents the cor- 
relation strength. The value of r^ • 100 is the percentage 
of variability in y associated with variability in x. 

Table 2 lists the Pearson linear regression correlation 
coefficients, r, for the 14 quantities explored. Columns in 
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Fig. 1. ZAHB model by Cassisi et al. (1999) fitted to NGC 1904 colour-magnitude diagram. Tlie effective temperature 
variation along the HB is also shown. 



Table 2 list the same quantities as in Table 1. In addition, 
the same correlation matrix has been computed for inter- 
mediate mctallicity clusters only (-1.8 < [Fc/H] < -1.3, 
Table 3). This allows to analyse the impact of the various 
correlations on a metallicity regime with a high sensitivity 
to any variation in the basic stellar parameters. 

We note that among the parameters studied here, only 
8 are measured independently, po is derived from /Uy(0), c 
and Tc', trh is derived from the My and r;j, etc. In general, 
correlations of any of the derived quantities with any of 
their constituent quantities or combinations do not pro- 
vide new information. On the other hand. Tables 2 and 3 
immediately suggest some interesting correlations that we 
will try to analyse next. 

— Metallicity: the first parameter 

The first correlation to be explored is the HB 

morphology-mctallicity dependence. As pointed out 
in the Introduction, metallicity is the so called first 
parameter regulating the extension of the horizontal 
branch, and its influence can be naturally derived from 
canonical stellar evolution models. Figure 3 shows the 
trend of log(TeffjirB) with [Fe/H]. Clearly, there is a 
correlation between both quantities in the sense that 
the less metallic the cluster is, the more extended its 
HB tends to be. However, the data indicate, as we 



already knew, that the variation of log(Teff//B) from 
cluster to cluster is not completely explained by the 
[Fe/H] parameter. This observational evidence is the 
core of the second parameter problem, mentioned in 
the Introduction. The value of the Pearson correlation 
coefficient for these two quantities is r ~ -0.54 (and r 
~ -0.57 in the intermediate metallicity regime), indi- 
cating that metallicity explains the ~ 30% of the total 
variation of log(TcffffB). A simple least-square fit is 
also plotted in Figure 3, giving a rms value of 0.19. 
The influence of the selected metallicity scale on the 
result has been analysed by repeating the calculations 
in the Carretta & Gratton (1997) metallicity scale. 
The value of the derived correlation coefficient be- 
tween [Fe/H] and log(Teffijs) slightly diminishes (r 
~ -0.47). 
— Total luminosity 

Probably one of the most interesting results of this sim- 
ple correlation approach is the finding of a clear corre- 
lation between HB morphology and total luminosity of 
the cluster. If no selection in metallicity is performed, 
this correlation is apparently slightly lower than the 
one observed with the first parameter. The Pearson 
coefficient relating these two variables is r ~ -0.48 and 
therefore, the variation of total luminosity would be re- 
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Table 1. Col. 1: cluster identification; Col. 2: logarithm of HB highest effective temperature; Col. 3: metallicity; Col. 
4: total cluster luminosity in V or absolute visual magnitude; Col. 5: coUisional parameter; Col. 6: logarithm of central 
luminosity density; Col. 7 central concentration; Col. 8, 9 and 10: distance from Galactic center, Galactic longitude 
and latitude (degrees); Col. 11 and 12: core radius and half-light radius in arcmin; Col. 13 and 14: logarithm of core 
relaxation time, and logarithm of relaxation time at the half-light radius, in logio (years); Col. 15: central surface 
brightness in V magnitudes per square arcsecond; Col. 16: Cluster relative age. 
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Table 2. Pearson linear regression correlation coefficients, r, for the 14 quantities in Table 1. 
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Fig. 2. Comparison between our Tes parameter, in its logarithmic form, 
(1993). 



and the Lt parameter from Fusi Pecci et al. 
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Table 3. Same as Table 2, but for intermediate metallicity clusters only. 



sponsible for the 23% of the HB morphology variation. 
However, if only intermediate metallicity clusters are 
considered, the correlation between log(TeffjifB) and 
My is as high as r ~ -0.77 (60 % of the total variation). 
Therefore, the influence of cluster total luminosity on 
the HB temperature extension seems more important 
than that of [Fe/H] at this metallicity regime! 
Graphically, this effect is shown in Figure 4, upper 
panel, where the correlation of log(TeffHB) with to- 



tal luminosity for 3 different metallicity intervals is 
presented. More luminous clusters, tend to have hot- 
ter (bluer) horizontal branches. The rms of a linear 
regression between both quantities is 0.20. On the 
other hand, if we calculate the Pearson coefficient be- 
tween total luminosity and log(TeffijB) only for the 
subsample of clusters with the more extended HBs, 
log(TeffijB) > 4.3, we find a much clearer correlation, 
reaching a value of r= -0.81. This means that two 
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thirds of the variation in the HB temperature extension 
can be explained by the variation of My for clusters 
with extreme blue horizontal branches. Bottom panel 
of Figure 4 shows again the trend of log(Teff//B) with 
metallicity, but using different point sizes depending 
on the value of the total luminosity. Clearly, some of 
the dispersion in the first parameter correlation can 
be explained by My. That is the case of the metal- rich 
blue tail clusters NGC6388 and NGC6441, which have 
My < -9, and are among the most luminous clusters 
in the sample. 

Luminosity is perhaps the most fundamental observed 
quantity characterizing a stellar system, and for a set 
of old stellar systems it is a good relative measure of 
its baryonic mass. Therefore, the observed trend seems 
to suggest that, for some reason, more massive clusters 
tend to have bluest horizontal branches. We will come 
back to this point in Section 7, where we will discuss 
its possible theoretical implications. 
On the other hand, correlation between HB morphol- 
ogy and total luminosity was previously noted by Fusi 
Pecci et al. (1993). However, they interpreted this re- 
sult as a consequence of the high HB morphology- 
cluster density correlation derived from their analysis. 
In this paper, we have chosen a different characteri- 
zation of cluster HB, the maximum temperature ex- 
tension, which shows a weak correlation with cluster 
density or even stellar collisions, as we will see later in 
this Section. 



Finally, we have explored the possibility that the 
Teff//B-My dependence could be a statistical effect, 
that is, the higher the number of stars in the cluster 
(and thus the more massive the cluster is) the higher 
the probability of finding hot HB stars. The so called 
"second parameter" could also be a mechanism with 
a low percentage of incidence, that would be only de- 
tected with high enough statistics. In principle, if this 
were true, we would expect the hot HB stars to be 
always a small percentage of the total number of HB 
stars, which is not always the case (see for example 
the cluster NGC6205). Nevertheless, in order to check 
the influence of statistics on the My parameter de- 
pendence, we have performed the following test: we 
have taken the photometry of one of the most mas- 
sive clusters in the sample, NGC 2808 (My = -9.36), 
whose extremely extended HB reaches a temperature 
of log(TeffijB)=4.568. The absolute total magnitude 
of the cluster was then reduced artificially by sub- 
tracting the corresponding percentage of stars from 
the photometry file, using a random selection proce- 
dure. Between the 34% and the 96% of the stars were 
removed for magnitude reductions between 0.5 mag 
and 3.5 mag. From the rc^sidting simulated CMDs, the 
highest temperature of the HB was measured follow- 
ing the same technique applied to the real clusters (see 
Section 3.1). This procedure was repeated 20 times 
for each simulated cluster magnitude. The results ob- 
tained are plotted in Figure 4, upper panel, where 
the points with error bars correspond to the mean 
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Fig. 4. Upper panel: correlation of log(Teff/fB) with total luminosity, for 3 different metallicity intervals: [Fe/H] < 
— 1.8 (open circles), —1.8 < [Fe/H] < —1.3 (filled circles) and [Fe/H] > —1.3 (filled triangles). Points with error 
bars are measurements from simulated clusters (e.g. the text). Bottom panel: HB morphology- metallicity correlation 
with different symbol sizes depending on total luminosity, My. Note the better correlation of clusters of intermediate 
luminosity. The smallest clusters ( My > -7.8 ) do not reach the bluest HB morphologies, and they all have log(Teff^fB) 
< 4.35. 
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log(TeffHB) value obtained for each cluster magnitude 
and its scatter. From our simidations, the HB temper- 
ature extension seems to decrease very little as cluster 
magnitude decreases (less than a 4% in all the mag- 
nitude range). On the contrary, the tendency for real 
clusters (points without error bars in Fig. 4), seems to 
indicate a steeper decrease of log(Teff/fB) with My- 
Therefore, the performed test proves that the role of 
statistics in the dependence of HB temperature exten- 
sion on cluster magnitude is very small. The origin of 
the log(TefF^fB)-My correlation could be a different 
physical cause, whose possible interpretation will be 
discussed later (cf. Section 7). 

— Collisional parameter 

Quite interesting is also the result for the collisional pa- 
rameter. The correlation between log(TeffHs) and the 
Tcoi is very low and probably inside the errors (r~0.14, 
about 2% of the total HB morphology variance). This 
value increases up to r~0.29, 8% of the total variance, 
if only intermediate metallicity clusters are considered. 
The relation between log(TeffjjB) and Fcoi is graph- 
ically presented in Figure 5, on the upper left panel. 
This result seems to suggest that even if the probabil- 
ity of stellar collisions is higher, the HB morphology 
is not affected. Close encounters and tidal stripping, 
suggested by Fusi Pecci et al. (1993) as a possible 
origin of bluer HB stars, do not seem to have a rel- 
evant role in HB morphology. Nevertheless, it is im- 
portant to note that the collision rate may have var- 
ied greatly through the clusters life time, especially 
in clusters which have undergone core collapse and re- 
expansion, and/or gravothermal oscillations. Although 
none particular correlation for core-collapse clusters in 
the sample has been noted, we should warn about the 
risk of overinterpreting the lack of a good correlation 
between log{TeSHB) and Tcoi- 

— Central density 

The scenario is very similar to that inferred from the 
log(TeffjyB) - Tcoj correlation. Figure 5, bottom left 
panel shows the trend of HB temperature extension 
with the central density of the cluster. Our analysis 
indicates a weak correlation between both quantities: 
r~0.16 for the complete clusters sample, and r~0.30, 
at intermediate metallicity [2.6% and 9.0% of the total 
log(TeffjjB) variance respectively]. Again, contrary to 
previous studies (Fusi Pecci et al. 1993), our sample 
indicates that cluster density is not a "good" second 
parameter, as it has a little influence on HB morphol- 
ogy. Nevertheless, an equivalent caveat to that of Tcoi 
applies here. The central density now may not be as 
relevant as the maximum density achieved in the past. 

— Other simple correlations 

Other quantities with which HB morphology seems to 
have a little but, maybe still significant correlations are 
the distance to the Galactic center, Rgc (rci0.19, 4% 
of the total logiTcSn b) variance) and half-light radius 
(r~0.18, 3%). The first one, may also be a secondary 



effect of the first parameter, as [Fe/H] has an already 
known trend with Rgc- 
- Age 

Finally, we have evaluated the influence of cluster age 
in a subsample of 47 clusters, in common with the De 
Angeli et al. (2005) data base (see Table 1 and Figure 
6). The relative ages of the 47 clusters go from 0.74 ± 
0.06 for the youngest cluster in the sample (NGC 362) 
to 1.085 ± 0.06 for the oldest one (NGC 7099). These 
relative ages are normalized to the average age of the 
most metal poor clusters {[Fe/H] < —1.7) as ex- 
plained by De Angeli et al. The derived Pearson cor- 
relation coefficient between log(TeffHB) and age is r= 
0.04 (0.2% of the HB morphology variance of the clus- 
ters subsample). This result confirms the fact, already 
pointed out by Rosenberg et al. (1999), that age can 
not be the only explanation for the second parameter 
problem. 

In general, cluster age most important dependences 
are those with metallicity (r= 0.38, 14%) and Galactic 
latitude (r= 0.42, 18%). 

However, it is interesting to point out that if we only 
consider intermediate metallicity clusters (filled circles 
in Fig. 6), the correlation between log(TeffjyB) and 
cluster age increases to r=0.76 (58%). This will be ex- 
plored later in more detail. 

5. The multivariate approach. 

While the simple approach of examining individual mono- 
variate correlations of log(TeffHs) with cluster parame- 
ters gives us a good first look at the problem of HB mor- 
phology, the complexity of the situation calls for a more 
sophisticated strategy. We are dealing with a multidimen- 
sional data set, in which sets of several observables may 
be connected in multivariate correlations. Simple, mono- 
variate correlations are only a special and rare case. In 
particular, the different trends of log(Teff//B) illustrated 
above indicate that the problem of GC horizontal branch 
morphology is intrinsically statistically multidimensional, 
and that must be attacked using a multivariate approach. 

First of all, we performed the Principal Component 
Analysis (PCA) on our data set, using all independent 
input variables. A code developed at the Instituto de 
Astrofisica de Canarias (lAC) by A. Aparcio has been 
used. We ignored the derived quantities, po, trc, ^rh, as 
they do not add to the dimensionality of the data man- 
ifold. In addition, among the three positional variables 
Rgc, I and b, we have selected only Rgc and b. In the 
same way, among the group of variables formed by Tcoi, 
fj,v and rc, we took only Tcoi, as it contains the other two 
parameters in its formula. The input data were renormal- 
ized by subtracting the mean and dividing by the sigma in 
each of the input variables. The number of the significant 
eigenvalues, i.e., those larger than expected from the mea- 
surement errors, gives the dimensionality of the data man- 
ifold. Each of the eigenvectors also accounts for a fraction 
of the total sample variance. It is generally agreed that 
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Galactocentric distance. Symbols indicate different [Fe/H] ranges, as in Figure 4. 



eigenvalues >1 are statistically significant, but somewhat 
lower ones may be as well, depending on how the data 
were normalized. 

Table 4 presents the eigenvalues (ei), fractional (Vi) 
and cumulative (d) contributions to the total sample 
variance, in percent, for the obtained PCA solution. 
Figure 7 shows the results of the PCA, as applied to the 
entire above set of independent input variables, in the form 
of correlation- vector diagrams. Usually, a steep drop in the 
successive eigenvalues or in the fractional contributions to 
the sample variance indicates where the number of statisti- 
cally significant dimensions stops, and where the noise be- 
gins. Nevertheless, the situation is not always so clearcut. 
In this data set, there could be at least four, but probably 
as many as six statistically significant dimensions or more. 

The first four eigenvectors account for the 78.8% of the 
total sample variance. They define a natural coordinate 
system for this data set. Projections of the input axes 
to the principal planes given by the eigenvectors (■^1,^2) 
and (Ci,'C3) are the correlation- vector diagrams as shown 
in Figure 7. In this representation, vectors corresponding 



to well-correlated variables define sharp or plane angles, 
on the contrary, uncorrelated quantities have orthogonal 
vectors. 

On the other hand, as in the previous Section, the sit- 
uation for intermediate metallicity clusters was also con- 
sidered, by performing PCA computations for the above 
described 8-parameters, including only clusters with -1.8 
< [Fe/H] < -1.3. In that case, the 100% of the data vari- 
ability is explained with only 7 eigenvalues (the first four 
eigenvalues account for the 82.8%), thus reflecting the di- 
mensionality decrease. 

Obviously, the situation is quite complex, reflecting 
the statistical multidimensionality of the entire manifold 
of globular cluster properties (e.g. Djorgovski & Meylan, 
1994). Given the limited data set, a more profitable ap- 
proach is to consider only a subset of variables. To illus- 
trate the point, we will consider only log(Teff//s), [Fe/H] 
and My. We find that for this data subset of three in- 
put variables (see Table 5), at least two statistical dimen- 
sions are necessary. The first two eigenvectors account for 
the 56% and the 34% of the total sample variance, re- 



Recio-Blanco et al.: HB Multivariate analysis 



11 



4.6 



4.4 



14.2 



0) 
O 



3.8 



n — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 



AO • 



o 



o 



O 



g g I I I I I I I I I I I I I I I I I I I I I I L_ 

■ 0.7 0.8 0.9 1 1.1 

Age 



1.2 



Fig. 6. Correlation of log(Tcff/fB) with globular cluster relative age. Symbols are the same as in Fig. 4 The derived 
Pearson correlation coefficient is r= 0.04 when all the clusters are considered, but it increases up to r=0.76 if only 
intermediate clusters (filled circles) are taken into account. 



i 


ei 


V. 


Ci 


1 


2.32 


29.0 


29.0 


2 


1.85 


23.1 


52.1 


3 


1.18 


14.8 


66.9 


4 


0.95 


11.9 


78.8 


5 


0.71 


8.9 


87.7 


6 


0.61 


7.6 


95.3 


7 


0.22 


2.7 


98.0 


8 


0.16 


2.0 


100.0 



i 


ei 


V, 


c. 


1 


1.68 


56.0 


56.0 


2 


1.01 


33.7 


89.7 


3 


0.31 


10.3 


100.0 



Table 5. PCA solutions for a subsample of 3 parameters 
: log(TeffaB), [Fe/H] and M^.. 



Table 4. Eigenvalues (ej), fractional (Vj) and cumulative 
(Cj) contributions to the total sample variance, in per- 
cent, for the obtained PCA solution for the set of 8 input 
independent parameters. 



spectively. The remaining 10% could be accounted by the 
errors, which are quite difficult to evaluate as explained 
in the previous sections. If this is true, a weighted vector 
sum of [Fe/H] and My vectors could correlate much bet- 
ter with log(Teffi/B) than [Fe/H] or My alone. This is 
shown in Figure 8 where the bivariate correlation involv- 
ing these 3 input quantities is: 

log{TeSHB) = -0.17-[Fe/H]-0.136- My + 2.84 

Clearly, the dispersion (rms=0.16) is now smaller than 

the one observed in any of the monovariate correlations 
for the total cluster sample. This bivariate correlation 
works specially well for clusters with log(TeffiifB) > 3.8 



(TeffffB > 6000 K), that is, clusters with at least some 
blue HB stars. In fact, if we perform the previous bivari- 
ate analysis considering only these clusters, we find a much 
clear correlation, with an rms=0.13, as shown in Figure 
9 ( we warn however of the smaller quantity of data in 
Fig 9). 

On the other hand, we have also explored the trivari- 

atc correlation of HB morphology with [Fe / H] , My and 
Age in order to possibly reduce the scatter in the previ- 
ous bivariate relation, due to possible effects of age. The 
analysis was performed for a subsample of the 47 clus- 
ters in common with De Angeli et al. (2005) which had 
[Fe/H] < —1.3. The corresponding PCA results are pre- 
sented in Table 6 and Figure 10. The third eigenvector 
significance, 14.8%, increases with respect to the combina- 
tion of [Fe/H] and My. Now, the first three eigenvectors 
account for the 93.0 % of the total variance. The resid- 
ual scatter is, therefore, of the order of the measurement 
errors. 
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6. Summary and discussion: theoretical As it has been extensively discussed by Djorgovski et 

implications. al- (1993) the GC manifold is rather complex, and this 

is perfectly true also when we focus our attention on the 
In this section, the results of the multivariate analysis pre- problem of the GC HB morphology and of its dependence 
sented above are summarized, trying to shed some light on the cluster parameters. In the previous section, we have 
in their possible theoretical meaning. shown that at least four parameters, but probably up to 
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Table 6. PCA solutions for the subsample of 4 parameters 
: log(TefrHs), [Fe/H], My and Age. 



eight, are needed to reproduce the HB extension in tem- 
perature. As Fusi Pecci, Buonnanno and collaborators (see 
for example Fusi Pccci ct al; 1993), among others, have 
been remarking in the last 20 years or so, there is no single 
"second parameter" which can explain the HB anomalies, 
but a combination of parameters. 

The present analysis, based on 54 GCs of the HST 
snapshot catalogue, has given a few interesting new results 
which are worth of further discussion. 

At least one new interesting conclusion can be derived: 
the importance of total cluster luminosity, and therefore of 
total m,ass, on the horizontal branch morphology. This ef- 
fect, combined with the first parameter, can explain prob- 
ably the major part of the Galactic globular cluster hor- 
izontal branch morphologies. More massive clusters (i.e. 
more luminous) tend to have more extended horizontal 
branches. To this scenario, we have to add the effect of 



age, evaluated here for 47 clusters in common with De 
Angeli et al. (2005). Lastly, the situation for intermediate- 
metallicity clusters has also been analysed, leading to a 
higher correlation between the temperature extension of 
the HB and My (60% of the total variation), with small 
increments of the Tcoi and the central density contribu- 
tions, that remain nevertheless inferior to 9%. 

One possible interpretation to the considerable influ- 
ence of My and therefore, of cluster total mass, on HB 
morphology can be derived from the paper by D'Antona 
et al. 2002. They analyze the consequences, on HB mor- 
phology, of helium variation due to self-pollution among 
globular cluster stars. Self-pollution had already been pro- 
posed as an explanation for the chemical inhomogeneities 
(spread in the abundances of CNO, O - Na and Mg anti- 
correlation) observed in GC members from the main se- 
quence to the RGB (see for example Gratton et al. 2001). 
The ejecta of massive asymptotic giant branch stars, which 
would be the origin of the self pollution, would not only 
be CNO processed, but also helium enriched. D'Antona et 
al. (2002) models take into account this possible helium 
enhancement with respect to the primordial value. They 
find that a spread in the helium content does not affect 
the morphology of the main sequence, turn off and RGB 
in an easily observable way. However, the difference in the 
evolving mass may play a role in the formation of blue 
tails, as higher helium stars would be able to populate 
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Fig. 9. Bivariate correlation of HB morphology with optimized combinations of [Fe/i?] and My for clusters with 
log(TefrHB) > 3.8. The plotted lines are the least square fit to the included data (solid) and to clusters with [Fe/iJ] < 
— 1.3 (dashed). Symbols are those of Fig. 4 



much bluer HB regions. If this is correct, self-pollution 
and so helium enrichment would be higher in more mas- 
sive clusters, as they would be able to retain the material 
from the ejecta better than less massive clusters. 

As already pointed out by Rosenberg, Recio-Blanco & 

Garcia-Marin (2004), M54, believed to be the remaining 
core of the Sagittarius dwarf galaxy could be another ex- 
ample of this scenario. More recently, the abundance anal- 
ysis of stars in the double main sequence of lo Centauri 
(Piotto et al., 2005) suggests the presence of two popula- 
tions of stars, one of which is strongly He enhanced. This 
could be another observational indication supporting the 
D'Antona ct al. theory in a very massive cluster. 

Finally, whatever is the theoretical interpretation of 
the data, a clear conclusion of this analysis is that the 
influence of My on HB extension seems to be as impor- 
tant as those of metallicity and age. Cluster total mass 
must be playing an important role in horizontal branch 
morphology. 
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